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Non-minimal supersymmetric models that predict a tree-level Higgs mass above the Minimal Su¬ 
persymmetric Standard Model (MSSM) bound are well motivated by naturalness considerations. 
Indirect constraints on the stop sector parameters of such models are signihcantly relaxed compared 
to the MSSM; in particular, both stops can have weak-scale masses. We revisit the stop-catalyzed 
electroweak baryogenesis (EWB) scenario in this context. We find that the LHC measurements of 
the Higgs boson production and decay rates already rule out the possibility of stop-catalyzed EWB. 
We also introduce a gauge-invariant analysis framework that may generalize to other scenarios in 
which interactions outside the gauge sector drive the electroweak phase transition. 


I. INTRODUCTION 

The origin of matter-antimatter asymmetry in the Uni¬ 
verse is a longstanding problem at the interface of particle 
and nuclear physics with cosmology. This issue cannot be 
addressed within the Standard Model (SM) and requires 
physics beyond the SM. One of the most interesting pos¬ 
sibilities, which has attracted much attention in recent 
years, is the electroweak baryogenesis (EWB) scenario 
(see [T] for review), where the baryon asymmetry is pro¬ 
duced during the electroweak phase transition (EWPT). 
This mechanism requires new physics beyond the Stan¬ 
dard Model (BSM) at the weak scale (^ 100 GeV) for 
two different reasons: 

• The SM with the Higgs boson mass above ^ 
80 GeV does not satisfy the Sakharov criterion [5] 
of departure from thermal equilibrium because the 
EWPT is a cross-over, rather than a strong 1st- 
order transition [3HS]. In order for EWB to be vi¬ 
able, one inevitably needs new bosonic fields that 
couple to the Higgs and significantly change the 
EWPT dynamics. 

• Although the SM does violate GP symmetry, the 
effects are highly suppressed at temperatures T ^ 
100 GeV by quark Yukawa couplings. Conse¬ 
quently, even if a SM universe admitted a strong 
Ist-order EWPT the produced baryon asymmetry 
is too small [US]. New sources of CP violation at 
the weak scale are required. 


Weak-scale supersymmetry (SUSY) is a well-motivated 
extension of the SM which can address both these 
problems. Electroweak gaugino/Higgsino [TUI [H] and 
scalar [T^ phases can provide the new sources of CP vio¬ 
lation (CPV), while the nature of EWPT can drastically 
change either in the presence of low-mass stops [nHH], 
or in models with extended Higgs sectors, e.g. additional 
singlet in the Next-to-Minimal Supersymmetric Standard 
Model (NMSSM) [TMl8] . In both cases, present and fu¬ 
ture experimental probes of BSM physics may provide 
conclusive tests. Present limits on the permanent electric 
dipole moments (EDMs) of the electron, neutron, and 
neutral atoms place stringent constraints on these CP- 
violating sources, and future EDM searches may probe 
the remaining CPV parameter space m- At the same 
time, searches for new scalar particles at the CERN Large 
Hadron Collider may uncover the ingredients needed for 
a first order EWPT. 

In this paper, we focus on the possibility that low-mass 
stops may give rise to the first order EWPT, the the so- 
called stop-catalyzed electroweak baryogenesis scenario. 
The basic idea is that a very light stop (m « 100 GeV) 
modifies the Higgs potential at finite temperatures via 
quantum loop effects, inducing a barrier along the Higgs 
direction between the electroweak symmetric and broken 
vacua and triggering a strongly Ist-order EWPT. It has 
been already shown in EOllll] that this scenario is no 
longer viable in the context of Minimal Supersymmetric 
Standard Model (MSSM), in light of the LHC measure¬ 
ments of the Higgs boson mass, production cross sections 
and decay rates (for pre-LHC works on the same subject 
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In this work we analyze the stop-catalyzed baryogen- 
esis in a more generic SUSY framework. Specifically, we 
consider adding small, hard-SUSY breaking terms to the 
Higgs potential. (Such terms can be dynamically gener¬ 
ated at the few-TeV energy scale; for examples, see [13- 
I26jl. The primary motivation for models of this type 
comes from considerations of naturalness [27]. The new 
potential terms can give tree-level contributions to the 
SM-like Higgs mass, allowing for a SUSY theory with an 
SM-like Higgs at 125 GeV independently of the stop sector 
parameters. Spectra with two relatively light stops can 
then still be viable, reducing hne-tuning. One may hope 
that the wide-open stop parameter space of these models 
may also allow them to accommodate the stop-catalyzed 
EWB scenario. This paper will explore whether this is 
indeed the case. 


Unfortunately, we find that even in this broader frame¬ 
work, there is still no parameter space which is compat¬ 
ible with both the strong 1st order EWPT and the LHC 
Higgs measurements. The key observation which leads 
us to this conclusion is that the light stop can change 
the order of the phase transition only if its coupling to 
the Higgs is close to its maximal possible value, which 
occurs when the light stop is nearly a pure gauge eigen¬ 
state (that is, U « or « in). On the other hand, 
in the small-mixing limit, loops of the light stop induce 
very large (~ 100%) shifts in the Higgs couplings to pho¬ 
tons and gluons, in contradiction with the LHC data that 
requires these couplings to be within 10 — 20 % of their 
SM values. It has been proposed in Ref. [5S] that this 
problem can be resolved if one assumes an appreciable 
invisible rate of the Higgs particle into neutralino pairs, 
which would compensate for the growth in hgg coupling 
due to the light stop. We show that this modification 
is also no longer compatible with the LHC Higgs data. 
Therefore, we conclude that stop-catalyzed EW baryoge- 
nesis is excluded by data over the entire parameter space, 
even if the MSSM Higgs mass constraint is removed. 


The paper is organized as follows. In Sec.jTTjwe briefly 
review the idea behind stop-catalyzed baryogenesis and 
describe in detail our framework of SUSY beyond MSSM. 
In Sec. |HI| we review the existing experimental con¬ 
straints on the light stops, including both direct LHC 
searches and indirect constraints from measurements of 
the Higgs properties. In Sec. W we discuss the condi¬ 
tions for a strong 1st order EWPT with an emphasis on 
gauge-invariance of the analysis. The framework that we 
introduce in this context may generalize to other scenar¬ 
ios in which new interactions outside the gauge sector 
drive the EWPT dynamics. Sec. jV] discusses some other 
technical details of our analysis, as well as its main re¬ 


sults. Finally, the conclusions are summarized in Sec. VI 


In the SM with the Higgs mass above > 80 GeV the 
EWPT is a cross-over W- Assuming that the phase 
transition occurs in a single step, any viable EWB sce¬ 
nario, thus, requires augmenting the Higgs sector with 
new bosonic degrees of freedom that couple strongly to 
the Higgs (for alternatives to the single-step phase tran¬ 
sition scenario, see, e.g. Refs. [2M35] b One generic pos¬ 
sibility is new scalar particles that generate a barrier be¬ 
tween the electroweak symmetric and broken vacua in 
the Higgs thermal potential via quantum loops, leading 
to a strongly 1st order EWPT. SUSY provides a natural 
candidate for such a particle: the stop. The coupling of 
the stop to the Higgs is predicted by SUSY and it is equal 
(up to the trilinear term) to the top Yukawa squared, an 
order-one number. It has been shown in unidiiEH] that 
light stops, with masses roughly in the 100... 120 GeV 
range, can trigger a strongly 1st order EWPT, if the stops 
are not heavily mixed and the tan/3 is large, > 10. The 
reason for this constraint is that mixing reduces the ef¬ 
fective coupling between the stop and the SM-like Higgs. 
Most of the work on this scenario assumed that the light 
stop is purely right-handed, to avoid introducing a very 
light sbottom. In this paper, we do not make this as¬ 
sumption, considering instead the most general stop sec¬ 
tor characterized by three parameters, the physical stop 
masses and , and the rotation angle from gauge 
to mass eigenbasis, Ot- (In some of the plots, we will find 
it useful to trade 9t for the mixing parameter Xt , defined 
as Xt = Mf^/mt-) Doing so will allow us to quantify 
the extent of mixing for which the strong Ist-order tran¬ 
sition is possible, as well as to contrast the parameter 
space required by the EWB scenario with that allowed 
by the LHC data. 

Discovery of the 125 GeV Higgs disfavors the pure 
MSSM, which predicts mh < rnz at tree level. The one- 
loop contribution to the Higgs mass-squared is given by 


Xnih = 


3yi 


—(r COS a mi log 

47r2 V 




ml 


( 1 ) 


where the ellipsis stand for the terms proportional to the 
Higgs mixing and higher-loop corrections [361138] . Very 
high stop masses are required to accommodate the mea¬ 
sured Higgs mass. For example, if « 100 GeV, as 
required in the stop-catalyzed EWB scenario, the sec¬ 
ond stop would need to have a mass of order 100 TeV 
to PeV, unless Xt > 1 TeV [31] • Such a spectrum would 
imply significant fine-tuning in the electroweak scale. In 
addition, the hierarchy of three or more orders of mag¬ 
nitude between the left-handed and right-handed squark 
soft masses would be challenging to explain from model¬ 
building point of view. Finally, even if these issues are 
ignored, the purely-MSSM version of the stop-catalyzed 
EWB is now in direct conflict with data, since the light 
stops loops give unacceptably large shifts to the hgg and 
/177 couplings MM- 
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In this paper we take a different approach, follow¬ 
ing Ref. [27] . We assume that additional TeV-scale 
physics beyond the MSSM provides a new tree-level 
contribution to the Higgs mass, reducing the need for 
large radiative corrections. There are many examples 
of such new physics, including (but not limited to) 
higher-dimensional F-terms |25[ I26j . or non-decoupling 
D-terms [531 If all non-MSSM states introduced 

by such models are assumed to be sufficiently heavy, 
with masses around I TeV or above, their effects at the 
^ 100 GeV scale, relevant for both the EWPT and Higgs 
LHC phenomenology, can be parametrized as new hard- 
SUSY breaking terms in the Higgs potential.^ These ef¬ 
fects can be parametrized by the generic two-Higgs dou¬ 
blet model (2HDM) potential: 

R = M^\Hu\^ + Mi\Hd\‘^ + {bHu ■ Hd + c.c.) + 

+ A 2 |H„P(i 7 „ • Hd + c.c.) + A 3 |i 7 „ni 7 dp + 

■ Hd + c.c.)^ -I- AsjiLu • Hd\^ -\- 
Xe\Hd\^{H^-Hd + c.c.) + ^\Hdf . (2) 

In the MSSM, the coefficients are given by 


As — — A 2 — A 4 — Ae — 0. (3) 

Additional BSM physics can change these relationships. 
While such corrections explicitly break SUSY, one should 
note that the potential Q is merely an effective descrip¬ 
tion valid at energy scales ~ 100 GeV, where SUSY is 
already broken. In the full UV theory, the corrections are 
generated by integrating out new physics at energy scales 
comparable to the soft SUSY-breaking scale, ^ 1 TeV. 

Not every term in the 2HDM potential can be easily 
generated by a UV-complete supersymmetric theory. For 
example, it difficult to see how any supersymmetric the¬ 
ory can generate non-zero coefficients A 2 , A 4 and Ag at the 
tree level. The other coefficients can be generated either 
by new F-terms or by new D-terms at multi-TeV scale. 
However, only dAi can give a substantial contribution to 
the SM-like Higgs mass without introducing new light 
states.^ For moderately large tan ,5, a tree-level 125 GeV 


^ Some SUSY models that lift the tree-level Higgs mass, e.g. 
nMSSM, introduce new states with masses ~ 100 GeV. These 
states can have an important effect on the EWPT 11711181 . We 
will not consider such models here. 

^ It has been shown in Ref. m that 5 A 5 alone can also solve the 
problem of the Higgs mass. However, this solution either works 
for tan/3 < 3, where there is no 1st order EWPT, or requires 
5 A 5 > 1, which can probably be UV-completed only with new 
light states. Therefore we will disregard this option. 


Higgs mass requires AAi < 0.1, which can be easily ac¬ 
commodated in the low-energy effective field theory and 
UV-completed at a scale of 2... 3 TeV [17]. Such hard- 
SUSY breaking term will be our basic assumption. Its 
presence completely removes the Higgs mass constraint 
on the stop masses and mixing. Apart from Ai, all other 
2HDM coefficients are assumed to be equal to the MSSM 
values. 

Parenthetically we note that UV-complete theories 
which induce a necessary dAi term will often also induce 
6 X 5 and AAy of the same order of magnitude. We will dis¬ 
regard these terms in our analysis because they have no 
significant effect on the SM-like Higgs mass. Moreover, 
if the 2HDM is not in the full decoupling regime, these 
terms usually only make the fit of the Higgs couplings 
worse. Therefore the limit where SX^j —>■ 0 can be re¬ 
garded as the best possible case for the light stop SUSY 
scenario. 

Two further simplifying assumptions will be made in 
our analysis. First, we assume that the 2HDM is in the 
decoupling regime, which for all practical purposes means 
^ 800 GeV. In this regime, the additional Higgs 
bosons beyond the SM-like 125 GeV state are sufficiently 
heavy to have no effect at the EWPT critical tempera¬ 
ture, Tc ^ 100 GeV.^ This assumption is strongly moti¬ 
vated by the agreement of the LHG Higgs measurements 
with the SM: any deviation from these limits would only 
further strengthen the LHC constraints on the scenario in 
the large tan/? regime of a supersymmetric theory. Sec¬ 
ond, we ignore the effect of all the superpartners, other 
than the stops, on the Higgs couplings and its thermal 
potential at T ~ Tc. This is motivated both by the non¬ 
discovery of superpartners at the LHC, and by the fact 
that the Higgs coupling to stops is the strongest among 
its couplings to the superpartners. Eor example, while 
electroweak gauginos and Higgsinos may be present at 
a few-hundred GeV scale (and indeed may be needed to 
provide CP-violating phases in the EWB scenario), their 
effect on the Higgs properties is subdominant to that of 
the stops. With these assumptions, both the EWPT and 
the collider constraints depend on just four unknown pa¬ 
rameters (m^^, rrij:^ , 0^, and tan /3) , enabling us to study 
the parameter space comprehensively and draw robust 
conclusions. 

The spectrum required for the stop-catalyzed EWB 
would necessarily require some degree of fine-tuning. At 
one loop, stop masses receive a quadratically divergent 
QCD radiative correction, cut off by the gluino mass rrig. 
Complete naturalness requires rrig ^ 2m^ in the case of 
Majorana gluino, or rrig ^ irrij: if the gluino is Dirac [42] . 
In either case, if the light stop is close to 100 GeV, gluinos 


® We assume in this paper m .4 Ri mjj Ri • While this is not 

true in generic 2HDM, and an appreciable hierarchy between 
these masses can have peculiar consequences on the EWPT (see 
e.g. |40lI41p . it is difficult to see how this hierarchy can be real¬ 
ized in a supersymmetric theory with no additional light states. 
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must occur below 400 GeV; this possibility is compre¬ 
hensively ruled out by the LHC data. Given the gluino 
bounds, a tuning of 0(10— 100) is required to accommo¬ 
date such a light stop. Our philosophy in this paper is 
to not be concerned about this; we would like to know 
whether or not the stop-catalyzed EWB is in agreement 
with the data, regardless of fine-tuning issues. 


III. CONSTRAINTS FROM DIRECT 

SEARCHES AND HIGGS MEASUREMENTS 

First, we consider direct searches for stops at the LHC, 
where the light stop in the 100 — 120 GeV mass range 
would be copiously produced. The sensitivity of the 
searches depends very strongly on the stop decay chan¬ 
nels. There are two possible options, corresponding to 
R-parity conserving and R-parity violating scenarios. In 
the R-parity conserving framework a light stop can decay 
into via an off-shell top (and possibly also an off- 

shell W). Alternatively, in the mass range where 3-body 
stop decay is prohibited, one can expect the decay mode 
i —>■ cx° to compete with the four-body decay. All these 
decay channels were studied by ATLAS and CMS collab¬ 
orations. The three- and four-body decays are searched 
for in monoleptonic channels [321 Ell • While most of the 
parameter space is excluded, the bounds are discontinu¬ 
ous, and become much weaker or even completely disap¬ 
pear near the borderlines between three- and four-body 
decays of the stop. Moreover, none of these searches is 
optimized for the four-body decays, and the constraints 
in this region are weak. This leaves a very light stop in 
the ~ 100 GeV mass range still a viable possibility. The 
searches for the stop decays in the two-body final 
state [ 321132 ] are more decisive and exclude the stops be¬ 
low the mass of ^ 200 GeV under the assumption that 
the branching ration (HR) into this mode is 1. However, 
the exact BRs in this part of parameter space depend on 
the neutralino mixing angles, stop mixture, and possible 
flavor violation in the scalar sector beyond MFV.'^ 

Another possibility is R-parity violation (RPV). In this 
case it is conceivable that the light stop is at the bottom 
of the spectrum and decays directly into SM states. Of 
course this possibility is excluded in the case of lepton- 
number violation. Searching for light stops in the baryon- 
number violating scenario, where the stop decays into 
two jets, with or without b-tag, is more challenging. 
However, recently GMS has excluded RPV stops below 
200 GeV [3Hjj rendering this option irrelevant for the 
stop-catalyzed baryogenesis. Even stronger constraints 
have been obtained from an ATLAS boosted search [32] , 
but they apply only for the b-tagged scenario. 


It has recently been claimed in Ref. dzl that the light stop sce¬ 
nario can be completely ruled out due to stoponium formation. 
It would be interesting to see whether a dedicated analysis by 
the experimental collaborations confirms this claim. 


To summarize, while direct constrains have already 
cornered the possibility of the light stops, the parame¬ 
ter space is not yet completely closed and the constraints 
are model-dependent. Therefore, we will now turn to 
analyzing indirect constraints, which are more robust. 

An important set of constraints on the light stop sce¬ 
nario comes from the electroweak precision measure¬ 
ments. Split scalar multiplets at the electroweak scale, 
such as stops and bottoms in the presence of mixing, con¬ 
tribute to the S and T parameters [50] . But the strongest 
constraints currently come from the measurement of the 
Higgs properties at the LHC. In particular, loops of very 
light stops significantly modify the coupling of the Higgs 
to the photons and gluons. If the additional Higgs bosons 
of the 2HDM are not too heavy, Higgs-fermion couplings 
can also be modified. In the decoupling approximation, 
which is almost always true in the SUSY context, and 
for moderately large tan/3 this effect dominantly modi¬ 
fies the Higgs couplings to the down-type quarks and the 
taus |2T] . 

Of course, the latter effect can be easily circumvented 
simply by decoupling the heavy Higgses. For example, if 
the heavy Higgs masses are around 800 GeV, we expect 
~ 3% correction to the h ^ bb rate, much too small to be 
detected with the currently available data. On the other 
hand, corrections to the couplings to photons and gluons 
are much harder to address. The stop loop contribution 
to the /i 77 and hgg couplings both scale in the small 
mixing regime approximately as 


9hgg/9 


SM 

hgg 


- 1 


1 f m]_ _ 

4 \ m? m? m? m? 

\ ti 12 ti 12 


(4) 


Without mixing, a 100 GeV stop produces an 0(100%) 
correction to the couplings, well beyond the 10 — 20% 
level allowed by current LHC data. The only way to 
cancel this correction is to turn on the mixing, which 
effectively suppresses the coupling of the light stop to 
the Higgs. But this is in conflict with the requirements of 
the EWB scenario, which requires a near-maximal stop- 
Higgs coupling and therefore small mixing. As we will 
see in Sec. jVj this tension cannot be reconciled with the 
current data. 

Before proceeding, let us note that the constraints from 
electroweak precision fits and the LHC Higgs measure¬ 
ments largely overlap, pointing to the same region in the 
stop parameter space [22H23] . In this region, called “fun¬ 
nel regime” or “blind spot” in the literature, the shift in 
the hgg/h'y^ couplings and the stop contribution to the 
T parameter are both minimized. (The stop contribution 
to the S parameter is small and plays a subdominant role 
in the fits.) This occurs at approximately 

2?tz+ 

sin(20) « - or Xt « m. -I- m. . (5) 

rrij:^ — rrij:^ ^ ^ 


In the limit relevant for the stop-catalyzed 

EWB, these conditions simply mean that the light stop 
coupling to the Hig^ vanishes, up to terms of order 
(see Eq. ® below). 
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IV. ELECTROWEAK PHASE TRANSITION 

The viability of EWB requires that the EWPT be 
sufficiently strong, that is, that the rate for baryon 
number changing sphaleron transitions inside the bro¬ 
ken phase, Esph, be slow enough to avoid washout of 
the baryon asymmetry. The rate Tsph is proportional 
to exp(—Egph/T), where Egph is the sphaleron energy at 
temperature T. The transition proceeds when T is below 
the bubble nucleation temperature, T^r, which is gener¬ 
ally just below the critical temperature Tc. The larger 
the magnitude of Egph/T/v, the more effective will be the 
baryon asymmetry preservation in the broken phase. 

In the context of perturbation theory, the computation 
of Tgph entails considerable conceptual and theoretical 
challenges [551156] . A particularly vexing one is mainte¬ 
nance of gauge-invariance. Loops containing gauge sec¬ 
tor degrees of freedom (gauge bosons, unphysical scalars, 
Fadeev-Popov ghosts) introduce gauge-dependence into 
the finite temperature effective action, Ses(T). Obtain¬ 
ing a gauge-invariant estimate of Psph at is possi¬ 
ble [55] , but doing so requires a level of care not typically 
followed in previous literature. 

Here we adopt a strategy that can be appropriate for 
the MSSM and other scenarios wherein gauge degrees of 
freedom play a subdominant role in generating the bar¬ 
rier between the symmetric and broken vacua. Specifi¬ 
cally, we truncate the one-loop effective potential I4ff(T) 
at second order in the electroweak gauge couplings g and 
g' while retaining while retaining terms to all orders in 
the top-quark Yukawa coupling, j/j. Doing so eliminates 
the gauge-dependence that first arises at 0 {g^) and that 
comes in tandem with the gauge-loop contribution to the 
barrier between the symmetric and broken-phase vacua. 
At the same time, it retains the gauge-invariant stop con¬ 
tributions to the barrier that enter first at 0 {yf) and that 
intuition tells us should dominate the phase transition 
dynamics. 

This intuition is based on the stop contribution to the 
daisy resummation term in Ves(T)-. 

APLsy(7^) = E , (6) 

i=l,2 

where Nc is the number of colors, “/i” generically de¬ 
notes the vacuum expectation values of the two neutral 
doublet Higgses, M^, (h) is the zero-temperature mass of 
stop eigenstate ti and M^, {h, T) is the corresponding fi¬ 
nite temperature mass. When the lighter eigenstate is 
essentially the right-handed stop, one has [13] 

(h, Tf D yfhl 1^1 - + ml + H,-^ (T), (7) 

where mg and mu are the left- and right-handed stop soft 
mass parameters, respectively, and Hj^ (T) is the one-loop 
thermal contribution to the stop mass-squared. Choos¬ 
ing mfj « —Hj^ (T) mitigates the screening of the stop 


contribution due to the daisy resummation. The result¬ 
ing approximate —Th^ term in the potential increases 
the barrier between broken and unbroken vacua, lowers 
Tc, and increases the ratio Egph/T as needed for baryon 
number preservation]!^. For T ^ 100 GeV, this choice 
leads to a lightest stop mass on the order of 100 GeV. 
Note that the coefficient of the stop-induced —Th^ term 
is enhanced by 2Ncyf- The gauge sector contributions, 
which are not included due to our truncation, carry no 
such enhancement. 

The requirements for effective baryon number preser¬ 
vation follow from solving the sphaleron equations of mo¬ 
tion and computing Egph- We observe that a consistent, 
non-trivial solution of these equations requires retain¬ 
ing gauge contributions to at least 0{g‘^) since the Higgs 
quartic self-couplings that enter the tree-level potential 
are as is the coupling between the gauge held and 

Higgs prohle functions [57]. In the present set-up, we for¬ 
mally retain all 0 {g‘^) contributions, but include none 
at higher order in g so as to maintain gauge invariance 
and consistency of the sphaleron equations of motion. In 
practice, for simplicity of numerical analysis we have not 
included the electroweak gauge boson contributions to 
the thermal masses that are also second order in g. We 
have estimated that doing so would result in shifts in the 
crucial quantity dehned in Eq. Q, by no more than 
10 %, leaving our conclusions unaffected. 

The resulting baryon number preservation criterion 
can be expressed as a condition on the ratio [58] 


where v{Tc) is the value of h{T) that minimizes 14ff (T) at 
the critical temperature. As discussed in Ref. [5S], there 
exist numerous sources of uncertainty in this condition, 
including the duration of the phase transition, the value 
of the baryon asymmetry at the start of the transition, 
the computation of the sphaleron fluctuation determi¬ 
nant, the origin of the unstable mode of the sphaleron, 
and neglected higher order loops [Ml EH [60]. Gonse- 
quently, the precise numerical results should be taken 
with a grain of salt. Nevertheless, we believe that our 
qualitative conclusions will not be altered, even taking 
into account significant uncertainties associated with the 
use of perturbation theory to analyze the phase transi¬ 
tion dynamics. To that end, we will show on our plots 
the contours of ^ = 0.5 and ^ = 0 (corresponding to 
the absence of a first order transition) to illustrate the 
potential impact of these uncertainties. 

V. RESULTS 

We performed a numerical scan over the four¬ 
dimensional parameter space outlined in sec.Jll] For each 
point in the scan, we evaluated consistency with the ex¬ 
perimental constraints by performing a ht to the LHG 
Higgs measurements, using the data set and the fitting 
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iH] = 100 GeV. tan /3 = 10 



400 600 


1000 1200 1400 1600 
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FIG. 1: Regions of stop parameter space allowed by the LHC 
Higgs measurements (blue - 67% CL and green - 95% CL) vs. 
the domain where the stop-catalyzed EWB can potentially 
be viable > 0) in pink. Mass of the light stop is fixed at 
100 GeV. The unphysical region (no solution for 9t) is shaded 
in purple. 


procedure of Ref. m, and the expressions for the cou¬ 
pling shifts from Ref. [52]. (For details of the fit, see Ap¬ 
pendix A.) We also evaluated the constraints from EW 
precision measurements; however, we find these to be 
consistently weaker that the Higgs fit constraints. Fur¬ 
thermore, for each point in the scan, we determined 
whether or not the EWPT is strongly Ist-order, using 
the procedure outlined in the previous section. 

We find that the Ist-order phase transition requires a 
very light stop, <110 GeV, independent of the other 
parameters. On the other hand, LEP-2 constraints imply 
> 100 GeV, confining this parameter to a narrow 
band. Within this band, the Higgs fit constraints on the 
remaining parameters vary only slightly with ^. In the 
plots below, we choose = 100 GeV as a representative 
value, but the picture that emerges from these plots is 
valid throughout the allowed range of m. Likewise, we 
fix tan /3 = 10 in the plots as a representative value. For 
larger tan (3 values our results stay almost independent of 
tan /3, while for lower tan /3 the EWPT becomes weaker 
while the Higgs constraints are largely unaffected. Again, 
the picture that emerges remains valid independent of 
tan/3. 

The main results of our analysis are summarized in 
Figs. Bi and The conclusion is clear: there is 
no overlap between the parameter space regions allowed 
by the Higgs fit, and those consistent with a Ist-order 
EWPT within a perturbative calculation. Thus, the stop- 
catalyzed EWB scenario is no longer viable. Perhaps 
the clearest way to understand this result is provided by 
Fig. HI Not surprisingly, the region of parameter space 


allowed by Higgs measurements is a band around the line 
Xt « rnf^: this follows directly from Eq. (|4]) in the limit 
. The crucial observation is that along the 
contours of constant EWPT strength Xt also scales 
linearly with mi^ . This is related to the fact that the ef¬ 
fective coupling between the lightest stop and the Higgs 
is given by 


-^e// ~ Vt 





(9) 


The thermal potential is determined almost exclusively 
by this effective coupling, so that constant-^ contours 
in the regime ^ correspond to a fixed ratio 
Xt/mi^. Grucially, a Ist-order transition is only possi¬ 
ble when the effective coupling is close to 1; specifically, 
Xt/mt^ 0.3 is required, as can be seen in Fig. This 
region does not overlap with the region ~ 1 al¬ 

lowed by the Higgs fits, regardless of the value of 
Incidentally, this argument provides a clear understand¬ 
ing of the results of Ref. |20| regarding the MSSM, where 
^ 100 TeV is required by the 125 GeV Higgs mass. 

Another useful representation of the same results is 
shown in Figs. an d§ where the Xt parameter has been 
traded for the stop mixing angle 9t. These plots make it 
clear for a 100 GeV light stop, the Higgs fits imply a tight 
relationship between 9t and which unfortunately is 
incompatible with the small-mixing regime required by 
the stop-catalyzed EWB. 



m- [GeV] 


FIG. 2: Constraints and regions with a Ist-order EWPT in 
the (mj^i^t) plane, for the light stop mass = 100 GeV. 
The region allowed by Higgs fits at 95% CL is shaded in blue. 
The blue line shows the contour of ^ = 0 and the red line 
shows the contour of ^ = 0.5. The region between the black 
contours is allowed at the 95% CL if a non-zero Higgs invisible 
width is included (etnv = 0.1 using the definitions of |61|L 
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tan B = 10, m- = 100 GeV 

h 



m-^ [GeV] 

tan B = 10, m- = 100 GeV 

h 



FIG. 3: Same as Fig. zooming in on the regions where the 
light stop is mostly II and Ir. 


Now we address an important potential caveat to the 
above argument. The Higgs coupling fit presented above 
is performed under the assumption of no new Higgs de¬ 
cay channels beyond those present in the SM. Lifting this 
assumption may be expected to relax the Higgs fit con¬ 
straints. In particular, the strong constraint on the gluon 
coupling comes from the agreement of the inferred pro¬ 
duction cross section in the dominant gg —)■ h channel 
with the SM. If the Higgs has a new decay channel, a 
production cross section in excess of the SM value can 
be accommodated without changing the observed event 
rates in any of the channels. In the stop-catalyzed EWB 
scenario, it is quite natural for the Higgs to have an ap¬ 
preciable decay width into light neutralino pairs, result¬ 
ing in an invisible final state. It was argued in Ref. [55] 
that this may revive the stop-catalyzed EWB scenario, 
even within the MSSM. 

However, since Ref. [55| appeared, dedicated experi¬ 
mental searches for Higgs invisible decays have been per¬ 
formed by both ATLAS and CMS, which strongly con¬ 
strain this possibility (55] [53] . Ref. m combined two 
these measurements (both of which had downward fluctu¬ 
ations) and estimated the allowed Higgs invisible branch¬ 
ing ratio as —0.18 ± 0.31, of course perfectly consistent 
with zero. The bound on /lyy coupling has also been im¬ 


proved significantly, so that the gluon coupling does not 
dominate the fit as strongly as in the early analyses of 
the LHC Higgs data. Due to both these factors, turning 
on the Higgs invisible width no longer improves the over¬ 
all Higgs couplings fit, and does not relax the constraints 
in the “right” direction for the stop-catalyzed EWB sce¬ 
nario. We explicitly illustrate this point in Figs. and [^ 
which shows that with the current data, allowing for the 
Higgs invisible width does not improve the fit. 


VI. CONCLUSIONS 

While stop-catalyzed EWB is a theoretically attrac¬ 
tive scenario, it has been known for some time that the 
LHC measurements of the Higgs mass and rates rule 
it out in the case of the MSSM. In this paper, we ex¬ 
tended the analysis to supersymmetric models with non- 
minimal contributions to the Higgs potential. Such con¬ 
tributions can lift the tree-level Higgs mass above the 
MSSM bound, thus eliminating one of the most stringent 
MSSM constraints on the stop sector parameters. One 
might hope that the newly expanded parameter space 
would include regions compatible with a stop-catalyzed 
EWB. We showed that this is unfortunately not the case. 
The reason for this is simple: the stop-catalyzed EWB 
requires a light stop ~ 100 GeV) with near-maximal 
coupling to the Higgs, since otherwise the effects of stop 
loops on the Higgs thermal potential are not large enough 
to trigger a strong Ist-order electroweak phase transi¬ 
tion. On the other hand, such a light stop is compatible 
with the LHC constraints on the hgg and /177 couplings 
only if its coupling to the Higgs is suppressed, and is 
far from maximal. We quantified these requirements and 
found no overlap between the parameter space regions 
with viable stop-catalyzed EWB and those with Higgs 
couplings compatible with the LHC data. The conclu¬ 
sion holds even in the presence of non-SM invisible Higgs 
decay channels such as /i —)■ 

An important limitation of our analysis is the assump¬ 
tion that no additional scalars beyond a single SM-like 
Higgs participate in the electroweak phase transition. 
This assumption can be violated in non-minimal super- 
symmetric models with extended Higgs sectors. For ex¬ 
ample, in the NMSSM-like scenario, an additional gauge- 
singlet scalar field can be active during the phase tran¬ 
sition, leading to novel possibilities such as a two-step 
phase transition. It has been shown that a strong Ist- 
order electroweak phase transition remains a viable pos¬ 
sibility in this class of models [T5lH8j . In these cases, 
the potential barrier necessary for the strong Ist-order 
transition typically arises from the tree-level multi-field 
scalar potential, rather than thermal loop effects, mak¬ 
ing them quite distinct from the stop-catalyzed scenario 
we focused on here. Still, these models demonstrate that 
even though the stop-catalyzed EWB no longer appears 
viable, other options for successful EWB exist in the su¬ 
persymmetric context. 











The analysis of this paper provides an illustration of 
the power of the Higgs data collected by the LHC to 
shed light on the electroweak phase transition, an im¬ 
portant event in the history of the Universe about which 
we currently have no direct information. The current 
10-20% accuracy of the Higgs coupling measurements 
is already sufficient to rule out one of the most popu¬ 
lar scenarios with a strong Ist-order phase transition. 
Many other models (both SUSY and non-SUSY) that al¬ 
low for a strong Ist-order phase transition predict more 
subtle deviations from the SM, and testing these models 
will require increasing the precision of the Higgs cou¬ 
pling measurements to the 1 % level and beyond, as well 
as measuring the Higgs cubic self-coupling [32l IMH 66 ] . 
The proposed electron-positron Higgs factories and a 100 
TeV proton-proton collider will be needed to comprehen¬ 
sively probe the viability of a strong Ist-order EWPT, 
and hence of electroweak baryogenesis. 
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Appendix A: Details of the Higgs Couplings Fit 

Constraints on the stop sector parameters were derived 
from a fit to the Higgs couplings measured at the LHC. In 
this Appendix we describe the fitting procedure in more 
detail. Since the procedure closely follows the prescrip¬ 
tion of Ref. m, one can view this Appendix as a short 
executive summary of this reference. 

The Higgs couplings taken into account in the fit are: 
hbb, hgg, h —>■ inv, hWW, hZZ, /177 and hrr. In the 
main part of our analysis, which assumes that the 2HDM 
is in the decoupling regime, the only deviations are in 
/177 and hgg and, in the case when Higgs decay to a neu- 
tralino pair is allowed, in h —>■ inv. The rest of the cou¬ 
plings only become important when we explicitly check 
the effects of 2HDM outside of the decoupling regime. 


which neccessarily triggers deviations in the couplings to 
the down type sector fermions and to the lesser extend 
to the gauge bosons. 

First, we slightly simplify the “proper” fit, by assuming 
that all error bars in the LHC measurements are Caus- 
sian (which does not lead to a significant loss of infor¬ 
mation). Second, we assume that all deviations from the 
SM-predicted values are small, namely = 1 -f £ 4 , where 
TiS are defined as 

t^higga = + 

V V V 

, (AI) 

and Ei <C 1. For the invisible rate, which is expected to 
vanish in the SM, we assume BR(/i —> inv) = emv To 
estimate the measured deviation of the 77 coupling from 
the SM we use the results of ATLAS [ST] and CMS [SB]. 

By combining these results we estimate e-y = 0.00 ± 0.14 
(of course we quote Icr uncertainties). The coupling hgg 
is extracted from all the exclusive (gluon fusion) decay 
modes, namely 77 [S7||SB], ZZ* [Sni[TD|, WW* [TTUT^] . 
and to the lesser extent from the fermionic channels. The 
average is eg = —0.13 ± 0.20. For the extraction of the 
bounds on the deviations of the other couplings, which 
are less important in our fit, the reader is again referred 
to Ref. [^. In the basis Ci = (et, eg, einv, £w, £z, ej, er), 
the correlation matrix between the different Ci’s, based 
on the theoretical calculation, is 



/ 1 

0.70 

0.04 

0.52 

0.38 

0.58 

0.59 \ 



0.70 

1 

0.43 

0.38 

0.11 

0.40 

0.52 



0.04 

0.43 

1 

0.46 

0.13 

0.40 

0.34 


p = 

0.52 

0.38 

0.46 

1 

0.44 

0.63 

0.45 

(A2) 


0.38 

0.11 

0.13 

0.44 

1 

0.42 

0.33 



0.58 

0.40 

0.40 

0.63 

0.42 

1 

0.54 



V 0.59 

0.52 

0.34 

0.45 

0.33 

0.54 

1 y 


For completeness we 

also 

quote 

here 

the allowed 

ranges 


of the other e’s, as calculated by m- 


et = -0.19 ±0.28, (A3) 

Gnv = -0.22 ± 0.20 , (A4) 

ew = -0.2 ±0.13, (A5) 

ez = 0.00 ±0.10, (A 6 ) 

= -0.03 ±0.17. (A7) 


Again, these channels were important in the fit only when 
we considered the non-decoupling limit of the 2HDM 
and/or the possibility of an invisible decay of the Higgs 
to light neutralinos. 
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